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Abstract. We present a derivation of the coherent backscattering spectrum 
from two two-level atoms using the diagrammatic approach, wherein the multiple 
scattering signal is deduced from single-atom responses, and provide a physical 
interpretation of the single-atom building blocks. 
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1. Introduction 

Coherent backscattering (CBS) of light emerges due to the constructive interference 
of multiply scattered counter-propagating waves surviving the disorder average 
pp. With optical waves, CBS was successfully observed using classical - e.g., 
polysterene particles [2] - and quantum - dilute cold atoms [3] - scatterers alike. 
A remarkable property of cold atoms is their nonlinear inelastic scattering induced by 
a powerful resonant laser field. Atomic saturation and inelastic scattering processes 
accompanying it were shown to affect the CBS interference from cold Sr (strontium) 
[4] and Rb (rubidium) [5] atoms, but a corresponding theory of CBS of intense laser 
light from cold atomic clouds is still lacking. 

The main challenge one has to deal with when describing CBS from cold saturated 
atoms can be briefly summarized as follows. A multiple scattering signal must be built 
on the basis of the accurately described responses of individual scatterers to a strong 
laser field. These responses can be found by solving the optical Bloch equations 
(OBE) [6]. However, a standard generalization of the OBE to the iV-atom case 
leads to a Lehmberg-type master equation governing the evolution of the reduced 
density operator of all atoms which are laser-driven and dipole-dipole interacting [7J|8], 
with the number of equations for the atomic averages growing exponentially with the 
number of scatterers. So far, such a master equation in the context of CBS has been 
solved only for N = 2 atoms [10]. 

Recently, we have initiated an alternative method of generalizing the OBE to the 
many-atom case which we call the diagrammatic approach to CBS [11] [T2l [13]. In 
its framework, the double scattering signal can be obtained from the solutions of the 
OBE for an atom subjected to a bichromatic classical driving. One component thereof 
represents the laser field and another one the field scattered from the second atom. 
Single-atom responses to a bichromatic field were evaluated non-perturbatively in the 
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Figure 1. Diagrammatic representation of the double scattering processes 
contributing to the background (a) and interference (b) spectra of CBS. The 
laser wave (thick arrows) of frequency ujl is scattered by the atoms (gray circles) 
into waves whose frequencies a/, oj" and ui^i may differ from u^. 



laser field and perturbatively - up to a second order - in the scattered field amplitude. 
Furthermore, by self-consistently combining single-atom responses (to which we will 
also refer as 'building blocks'), we were able to derive analytical expressions for the 
background and interference spectra of double scattering. The spectra thus evaluated 
were rigorously shown to be equivalent to that deduced on the basis of the two-atom 
master equation [Tl] . 

The motivation of the present contribution is twofold. First, previously we focused 
on a discussion of the inelastic building blocks and their self-consistent combination 
into double scattering diagrams [THEI]. These blocks determine the CBS signal only 
in the case of a very strong laser driving. In the present contribution we will consider 
a general case and present a diagrammatic derivation of the CBS spectra for arbitrary 
intensity of the laser field. Second, the single-atom building blocks have not been 
given a physical interpretation. Here, we furnish the single atom building blocks with 
a physical interpretation and establish a close connection with the results of Mollow 

EM- 

The paper is organized as follows. In the next section we present the building 
blocks contributing to the double scattering background and interference spectra. 
Thereafter it is shown how these blocks can be evaluated by solving the OBE under 
bichromatic driving. Thereby we establish the connection to the method used by 
Mollow in [T5i . In Sec. [3] we formulate rules for combining single-atom building blocks 
into double scattering diagrams, present the full set of diagrams contributing to the 
elastic and inelastic spectra, and give the explicit expressions thereof. We conclude 
our work in Sec. |U 

2. Single-atom building blocks 

2.1. Graphical representation of the elastic and inelastic building blocks 

To be self-contained, we will briefly outline here the main idea of the diagrammatic 
approach to CBS of laser light from two two-level atoms which was presented in 
detail in [T3] . The fundamental double scattering processes surviving the disorder 
averaging and contributing to the CBS background and interference signals are shown 
in Fig. [TJ We assume that the laser field may be sufficiently strong to saturate the 
atomic transitions. Accordingly, the frequencies a/, w" and wj) of the waves scattered 
by the atoms towards each other and a detector can, but need not, differ from the laser 
frequency w^, i.e., correspond to inelastic scattering. The co-propagating positive and 
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Figure 2. Decomposition of scattering processes at each of the atoms (gray 
circles) from Fig.[T]as a sum of elastic (open circles) and inelastic (hatched squares) 
building blocks. Processes contributing to (a) background, and (b) interference 
intensities. Frequencies of the probe and scattered fields will be defined after the 
expressions for the building blocks on the right hand side, as well as the rules for 
combining the building blocks are known. 



negative frequency amplitudes (solid and dashed arrows, respectively) contribute to 
the background double scattering intensity, see Fig. QJa), which is independent from 
the observation direction, whereas the counter-propagating ones contribute to the CBS 
interference, see Fig. HJb) ■ We will be interested in finding the frequency distributions 
of the background and interference spectra (intensity vs. wu) as a function of the 
driving field parameters (such as the Rabi frequency and the offset from the atomic 
transition frequency) . 

In the framework of the diagrammatic approach [TTJ [JH [T3], these frequency 
distributions can be derived on the basis of single-atom building blocks. The latter 
represent spectral responses of an atom subjected to a classical bichromatic driving 
field. The first component thereof corresponds to a laser field of arbitrary strength, 
whereas a second, weak, component describes the far-field scattered by the other 
atom. Following the nomenclature used in laser spectroscopy [51 116). we will refer to 
the laser and weak field components as the pump and probe fields, respectively. While 
the classical description of the laser field is common, applying the same description 
to the atomic radiation exhibiting photon antibunching [17] is, in general, wrong. 
However, in the dilute regime, when the double scattering originates from exchange 
of a single photon, the validity of the semiclassical treatment of the atom-probe field 
interaction does not contradict the nonclassical character of the scattered field PJ2[ Q/3] , 
which manifests itself only in correlations between at least two photons. Furthermore, 
the validity of the classical ansatz for the probe field was proven analytically by 
establishing the equivalence between the results of the diagrammatic and master 
equation calculations [T%| . 

The classical ansatz for the fields exchanged between the atoms allows us to 
consider two atoms with their incoming and outgoing fields separately from each other. 
The left hand sides in Fig.^a) and (b) show the decomposition of the background and 
interference contributions to CBS from Fig.[]Ta) and (b), respectively, into single-atom 
blocks. In order not to overburden diagrams, in Fig. [5] only fields scattered by the 
atoms are depicted, but here and henceforth one should remember that the atoms are 
also laser-driven and, furthermore, the effect of the laser field on atoms is accounted 
for non-perturbatively. Note that the arrows are not labeled by their frequencies in 
Fig. [2J The latter will be defined in the course of the subsequent analysis. 

In general, for any finite value of the saturation parameter, there are non-zero 
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cross-sections for elastic and inelastic scattering of photons by a single laser-driven 
atom [BJ. The same also holds in presence of additional probe field(s). In fact, one 
of the purposes of the diagrammatic approach is to calculate elastic and inelastic 
responses of a single laser-driven atom subject to probe fields. 

Accordingly, the right hand sides of Fig. [2] represent decompositions of the 
total single-atom responses into elastic (blank circles) and inelastic (hatched squares) 
building blocks. Different shapes of the blocks (i.e., circles or squares) emphasize that 
the corresponding expressions result from different equations of motion (for the atomic 
dipole averages and temporal correlation functions, see |Appcndix A[ ) describing the 
elastic and inelastic scattering, respectively. Different colours refer to photons emitted 
at the laser frequency (blank), or at the frequency different from the laser frequency 
(hatched), as will be explained in Sec. 12.2.21 

There is an important difference between the diagrams in Figs. Q] and [2J While 
the former diagrams depict general (background and interference) double scattering 
processes for a random configuration of atoms, the latter ones show single-atom blocks 
which, by construction, only produce those background and interference contributions 
which automatically survive the disorder average. In accordance with this, all arrows 
in Fig. [2] are for convenience oriented along the horizontal line. 

Prior to presenting the explicit expressions for the building blocks on the right 
hand side of Fig. [5J we will show that the elastic blocks can be decomposed further. 

2.1.1. Elementary 'elastic' building blocks Each of the building blocks in Fig. [2] 
contains two outgoing arrows. We will now express all elastic blocks (circles) as 
combinations of the blocks which each contain only one outgoing arrow. We will 
refer to such building blocks with one outgoing arrow as 'elementary' ones. 

In the case when there are no incoming arrows, see Fig.^a), the outgoing arrows 
correspond to the elastic intensity of light scattered by the laser-driven atom. It is well- 
known [BJ that the elastic intensity is proportional to {r7 + er~)(°' cl ) = (cr + )' ' (er~)' ', 
that is, given by the product of the expectation values of the atomic dipole operators 
for an atom driven only by the laser field (we will supply the expectation values in 
this case by the superscript (0)). We represent this equation graphically in Fig.[3][a). 
There, a circle with the outgoing solid (dashed) arrow corresponds to (cr^)^ ((a + )^), 
respectively, whereas the product of these two amplitudes is depicted by symbol 'x'. 

It is easy to generalize this result to the case when there are incoming probe-field 
amplitudes. In this case, the elastic building blocks can be expanded into sums of 
products of two elementary amplitudes, with the number of terms in this expansion 
being equal to the number of ways in which the incoming arrows can be distributed 
among the two circles (see Fig. E^b-d)). 

2.1.2. 'Inelastic' building blocks Concerning the inelastic building blocks of Fig. [2] 
(hatched squares), they result from the non-factorizable part of the atomic response, 
i.e., (a + a~) — (o- + )(a~), see equation (J6J) below. Unlike the elastic blocks, these 
blocks, which exhibit two outgoing arrows corresponding to a + and cr~, respectively, 
cannot be factorized into a product of blocks with only one outgoing arrow. 

2.2. Calculating building blocks 

In Sec. 12.11 we presented the disorder averaged elastic and inelastic single-atom 
building blocks which contribute to the double scattering background and interference 
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Figure 3. Decomposition of the elastic scattering processes into elementary 
processes, (a) the elastic intensity of light scattered by a two-level atom is equal 
to (al): the product of the amplitudes for the elastic scattering of positive and 
negative frequency amplitudes described by circles with outgoing solid and dashed 
arrows; (b) in case of two incoming probe fields there are four ways, (bl)-(b4), to 
distribute the probe fields among the irreducible blocks; (c), (d) decomposition of 
the elastic responses in case of one incoming probe field into (cl), (c2) and (dl), 
(d2), two elementary blocks, respectively. 



spectra of CBS. All these blocks represent a single, two-level, laser-driven atom which 
additionally receives none, one or two classical probe-field amplitudes. We will next 
show how to evaluate these blocks using the OBE under bichromatic driving. 

2.2.1. Single-atom optical Bloch equations under bichromatic driving A classical 
bichromatic driving held reads 

£{t) = Sle luJLt + £ L e- luJLt + e*e iu ^ + ee~ l ^\ (f ) 

where £l, e are complex amplitudes, and lul, uj p the frequencies of the pump 
and probe held, respectively. It is easy to show that the optical Bloch vector 
(a) = «er->> + » z », where a~ = |0)(1|, a+ = |1)(0|, a z = |1)(1| - |0)(0|, and 
|0), |1) the ground and excited states of the atom, obeys the following equation of 
motion written in the frame rotating at the laser frequency: 

(c?) = M{a) +L + v e lut /S. {+ \a) +v*e- lu}t ^-\a), (2) 

where M denotes the optical Bloch matrix for an atom driven by the laser field: 

/ -j + iS -iVt/2 \ 
M = -7 - iS ifl*/2 , (3) 

\ -m* m -27 / 



with 6 = ujl — wo the laser detuning from the atomic transition frequency ujq, and 7 
half the spontaneous decay rate, n — 2£l/H and v = 2ed/h are the Rabi frequencies 
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of the pump and probe fields, respectively, with d the (real) matrix element of the 
atomic dipole transition, uj = uj p — uj^ is the detuning between the probe and pump 
frequencies. Finally, the matrices 

/ -i/2 \ / \ 

A(") = , A( +) = i/2 \ (4) 
\ i 0/ \ -i J 

describe the coupling of the atom to the positive and negative frequency components 
of the scattered field, respectively, and 

L = (0,0,-2 7 ) T (5) 

is a constant vector. 

2.2.2. Elastic building blocks It is straightforward to establish the correspondence 
between the elementary elastic building blocks of Fig. [3] and solutions of equation 
We will start by noting that if v = v* = in equation ([2]), it reduces to the standard 
OBE having the steady-state solution (o^ -*. For non-zero probe fields, there appear 
corrections to the unperturbed components of the Bloch vector as well as oscillations 
at the frequencies ±w (see | Appendix A"| . Depending on whether a single incoming 
probe field represents a positive- or negative- frequency amplitude at frequency uj, the 
corrections will be denoted (<r(w)}(~) and (<r(w))( + ), respectively. A second-order 
correction, which corresponds to one solid and one dashed incoming arrow at the 
frequency u>, will be denoted as (ct(cj)}( 2 ). Finally, we recall that the emission of the 
positive-(negative-) frequency amplitudes is associated with the expectation values of 
the atomic lowering (raising) operators [6]. Summarizing the above, we express the 
elementary elastic building blocks in terms of the solutions of equation @ as shown in 
Fig. [4] From now on we use the convention not to label the incoming (outgoing) waves 
which are exactly on resonance with the laser frequency. Furthermore, we will denote 
elastic building blocks with outgoing wave at the laser frequency by a blank circle, 
whereas a hatched circle corresponds to an outgoing wave at a frequency uj ^ u>l (as 
shown below, the latter requires the presence of an incoming wave with frequency uj 
or — uj). 

Diagrams (a) and (b) of Fig. [4] describe positive and negative frequency 
amplitudes, respectively, of the elastic scattering by a laser-driven atom. Clearly, 
their product is nothing but the diagram (al) in Fig. [3] Diagrams (c)-(h) include 
one or two incoming arrows at the frequency uj. Atomic response functions depicted 
by these diagrams were studied by Mollow [TS] in the context of the weak probe 
absorption or amplification by a laser-pumped two-level atom (see | Appendix B"j ). As 
follows from the analysis presented in |Appcndix B[ the processes (c) and (d) are 
correspond to the amplitudes of stimulated emission at the probe field frequency, 
and hence, their outgoing arrows are labeled by the same frequency values as the 
incoming ones. In accord with the convention adopted above, these building blocks 
are hatched, since the frequency of the incoming wave - and therefore also of the 
elastically scattered outgoing wave - differs from the laser frequency. Further on, 
the processes (e) and (f) correspond to a nonlinear phase conjugation whereupon 
the positive- (negative-) frequency amplitude at the frequency uj is scattered into a 
negative- (positive-)frequency, that is, conjugated, amplitude at the frequency —uj, via 
a nonlinear mixing with the laser wave described by the nonlinear atomic susceptibility. 



Diagrammatic approach to coherent backs cattering of laser light by cold atoms: Double scattering revisited! 

(a) (b) 

= <<r><™ 0---> = <<7 + > <0) 

(c) (d) 

A.^-^ = <a-(w)}<-> -%#-% = (cr + (oo)> <+) 



00 ^ 



(e) (f) 

(+) = /rr+^\(-) 



>®^» =(ct-(uu)> <+) A#---> = <a + (uu)> 

(g) (h) 

={cj-(o))) (2) ftO > = <tf + ("»> (2) 

Figure 4. Graphical representation of the elementary diagrams for elastic 
scattering, and their mathematical expressions in terms of solutions of equation J2} 
for the optical Bloch vector. Waves at the laser frequency are not labelled. 
Left and right columns describe positive- and negative-frequency amplitudes, 
respectively, of: (a) and (b) - elastic scattering in absence of the probe field; 
(c) and (d) - amplitudes of stimulated emission at the probe field frequency ui by 
a laser-driven atom; (e) and (f) - phase conjugation via nonlinear wave mixing 
of the pump and probe fields; (g) and (h) - amplitudes of stimulated emission at 
the laser frequency by an atom subjected to a positive- and negative-frequency 
probe field at the frequency uj. 



Using the nonlinear optics terminology, one can regard [a~ ~{lS))(~\ (<t + (and 
their complex conjugate quantities) as, respectively, the effective (that is, dependent 
on the laser pump) linear and nonlinear susceptibilities of a two-level system with 
respect to the probe field [IS] ■ We note that the relevance of nonlinear susceptibilities 
for CBS from saturated atoms was mentioned before [19] ■ Finally, the processes (g) 
and (h) describe stimulated emission of the amplitudes at the pump field frequency in 
the presence of the probe field. 

We conclude this part by a short remark. The diagrams of Fig. |4] describe all 
possible elementary building blocks that appear in Fig. |3] Consequently, solutions of 
equation @ for the Bloch vector to second order in the probe-field amplitude allow 
us to evaluate completely the 'elastic' responses on the right hand side of Fig. [2J 

2.2.3. Inelastic building blocks Let us now address the inelastic building blocks. 
They are deduced from the solutions of equation @ for the atomic dipole temporal 
correlation function [T31 [T^] (see also |Appcndix A[ ) , 

g(t u t 2 ) = (<7+(t 1 )<7-(t 2 )) - (<7+(tx))(a-(i 2 )), (6) 

in the steady-state limit when g{t\,t2) — g{t\ — £2), and subsequent Laplace 
transformation to the frequency domain. 

The Laplace transformed solutions of equation ©, with the probe field amplitudes 
v and v* set equal to zero, together with the linear (oc v,v*) and quadratic (oc 
|v| 2 ) corrections, yield the inelastic building blocks shown in Fig. [SJ The physical 
interpretation of the functions P'°'(u;), P^^w^v) and P( 2 \us,v) is analogous to the 
interpretation of the building blocks of Fig. 01 and merely generalizes it to include the 
multi-photon scattering of the laser photons. 
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Figure 5. Graphical representation of the inelastic building blocks and their 
mathematical expressions (see equation ||A,4|| ). (a) P^(uj) - the Mollow triplet; 
(b) P( 2 \uj,v) - correction to the Mollow triplet due to the probe field at the 
frequency ui; (c) v) and (d) P(~ )(o), u) — inelastic wave mixing of the 

pump and probe fields resulting in scattering of the waves neither of which is in 
resonance with the laser or probe waves. 

Diagram (a) in Fig.[5]represents the inelastic spectrum of a two- level atom subject 
to monochromatic driving. That is, the function P^{uj) is the resonance fluorescence 
spectrum known, for f2 ^ 7, as the Mollow triplet [2D]. Diagram (b) represents the 
correction to the Mollow triplet p( 2 )(w, v) proportional to the intensity of the probe 
field at the frequency lj. This diagram can be regarded as a generalization of diagrams 
(g) and (h) of Fig.|4j which accounts for the inelastic scattering. Finally, the functions 
P^(u, v) (see Fig. [5jc,d)) represent inelastic wave mixing processes resulting in the 
emission of two amplitudes which are off-resonant with either the laser or the probe 
fields. These diagrams should be contrasted with the diagrams (c,d) in Fig. [3] Indeed, 
in the latter diagrams, the two outgoing arrows are on-resonant with the probe and 
laser fields. 



3. Combining building blocks 

3.1. Rules 

In the previous Section we defined all single-atom building blocks. We recall that 
the elastic and inelastic building blocks from which the background and interference 
contributions are constructed are shown on the right hand sides of Fig.[2jja) and (b), 
respectively. The four elastic blocks in Fig. [2] can be further decomposed according to 
Fig. [3] with the elementary diagrams defined in Fig. [4] whereas the inelastic blocks and 
their corresponding spectral response functions are shown in Fig. [5] So, all diagrams on 
the right hand side of Fig. [2] are defined, and we need to introduce rules for combining 
them. 

These rules can be formulated as follows: 

• The ladder and crossed contributions are obtained by combining the building 
blocks of atom 1 with that of atom 2, on the right hand side of Fig. [2ja) and 
(b); each combined two- atom diagram is proportional to the factor g 2 , where 
g = (k^i)^ 1 scales as the far-field dipole-dipole coupling between the two atoms 
separated by the mean distance £. 

• The elastic building blocks for atoms 1 and 2 are expanded according to Fig. [3] 

• Each outgoing solid (dashed) arrow of atom 1 is merged with each incoming solid 
(dashed) arrow of atom 2. As a result, one obtains sums of diagrams containing 
two intermediate and two outgoing arrows. The intermediate arrows' frequencies 
are defined self-consistently: the frequencies of the incoming and outgoing arrows 
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Figure 6. Background contributions to double scattering, (a) Fully elastic 
contribution originating from both atoms scattering elastically; (b) inelastic 
contribution due to second atom scattering inelastically; (c) mixed contribution 
including diagrams (cl), (c2) corresponding to elastic scattering, and (c3), (c4) 
corresponding to inelastic scattering; (d) fully inelastic contribution. 



must coincide and, moreover, respect the form of the response functions presented 
in Figs. 21 [5j As a result, for example, there is only one independent inelastic 
frequency describing intermediate amplitudes. 

• The response functions associated with the single atom building blocks in 
the resulting double scattering diagrams must be multiplied. The symbolic 
expressions for these functions are given in Figs. [4] and [5l 

• If an intermediate arrow representing the inelastic probe field (say, at frequency 
w) changes its frequency then the expression corresponding to that diagram is 
integrated over this frequency as /_ du>/(2ir) (note that u) — u p — u>l, and can 
take on negative values). 

• Diagrams containing blocks connected by intermediate counter-propagating 
arrows with no outgoing arrows must be excluded. Such diagrams correspond to 
a photon cycling between the atoms without contributing to the detected signal. 
Exactly the same rule exists in the nonlinear diagrammatic theory of CBS in a 
Kerr medium [21] . 

Using the above rules, we can represent the background and the interference 
contributions as sums of the diagrams shown in Figs. [6] and [7J respectively. When 
writing down the explicit expressions for different contributions, we will omit the 
overall factor g 2 . 



3.2. Background contributions 

3.2.1. Elastic spectrum By definition, the elastic background spectrum L c \(v) = 
L e \5{v) exhibits a single 5-peak precisely at the laser frequency, corresponding to 
v = in the rotating frame. Hence, it is given by the sum of all diagrams with 
unlabeled signal (outgoing) lines in Fig. HO that is, diagrams (a), (cl), and (c2). Note 
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that diagram (a) is to be unfolded into a sum of products of the following diagrams 
of Fig. [3J (al)(bl), (al)(b2), (al)(b3), and (al)(b4). Taking the expressions for the 
building blocks from Figs. 0] [SJ and applying the above rules for combining these 
building blocks, the elastic background intensity reads: 

+ <a + }( >(O (0 V + } (t V~(0)) (2) 
+ <a+)( )(O(°)( ( 7+(0))(+)<a-(0))(-) 

+ (a+)( )(O(°)(a-(0)}M(a+(0)) M (7a) 
+ /"JpMN^M)* 2 '^)' ) (7b) 

+ /^^ (0) H(^H) (2) {- + ) (0) , (7c) 

where equations (ITal) . f|Tb[) and (ITcl) correspond to diagrams (a), (cl) and (c2) in 
Fig. |51 respectively (see |Appendix A| for the explicit evaluation of the right hand side 
of equation (J7J)). 

3.2.2. Inelastic spectrum As evident from Fig. [6j its remaining diagrams describe 
the inelastic spectrum 

L- mel {v) = ( ff + ) (0) (a-)(°)Fl 2 )(0,,) (8a) 
+ P {0) {v)(o--{v)) { - ) {<t + {v)) (+) (8b) 
+ P(°) {v) (<t+ (i/)) (-) (a- (i/)) (+) (8c) 

+ /°° ^ (81 HP (2) (^), (8d) 

where equations (|8a|). (|8b|) . (|8c|) . and (j8d| correspond to diagrams (b), (c3), (c4) and 
(d) of Fig. ITU and the (signal) frequency uj was for convenience relabeled in lines 2 and 
3 of equation © to v. 

3.3. Interference contributions 

3.3.1. Elastic spectrum By the same argument as in the case of the background 
diagrams, we identify the diagrams (a), (bl), and (cl) of Fig. [7] as contributing to 
the elastic spectrum. In the present case, the fully elastic contribution represented by 
diagram (a) is explicitly decomposed into diagrams (al), (a2), and (a3) corresponding 
to a combination of diagrams (cl)(d2), (c2)(dl), and (c2)(d2) of Fig. |3l respectively. 
The product of diagrams Fig.|3jcl) and (dl) does not contribute to Fig. [7](a) because 
it corresponds to a photon cycling between the atoms (see Sec. 13. 1| . 

Note that the squares which are hatched in Fig. EJb), (c) become blank in 
diagrams (bl) and (cl). This is precisely the consequence of the self-consistent 
combination of diagrams: the connection of the phase conjugation diagram Fig. |U[e), 
for example, to Fig. EJd) enforces the occurrence of two frequencies uj and — w, and 
therefore v — 0, in the building block Fig. [5fd) . Thus, a blank square denotes an 
inelastic process, where the frequency of one of the outgoing arrows equals the laser 
frequency. 
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Figure 7. Interference contributions to double scattering, (a) Fully elastic 
contribution; (b) and (c) mixed contributions including diagrams (bl), (cl) 
corresponding to elastic scattering, and (b2), (c2) corresponding to inelastic 
scattering; (d) fully inelastic contribution. 



After these remarks, we present the explicit expression for the elastic interference 
spectrum: 

C cl = ( ( 7+>(°){ct+}( )(ct-(0))(+){ ( 7-(0))(-) (9a) 
+ (O (0) (O ( °V + (0)> ( ~V + (0)> (+) (9b) 
+ (a-}(°V + ) (0 V~(0)> ( ~V + (0)> (+) (9c) 

+ r ?(^ + ) (o) (a-(- W ))(+)p(-)( w ,0) (9d) 
J-oo 2?r 

^( ff -)(0)( (J +(- w ))(-)p(+)(a; ) 0), (9e) 

where the subsequent lines of equation © correspond to diagrams (al), (a2), (a3), 
(bl), and (cl) of Fig. [71 respectively. 

3.3.2. Inelastic spectrum The three contributions to the inelastic spectrum are 
depicted by diagrams (b2), (c2) and (d) in Fig. [7J In diagram (d), both the solid and 
dashed lines corresponding to intermediate amplitudes that change their frequencies. 
To write down the expression for this diagram, we note that the outgoing inelastic 
photon frequency is v. Since the frequencies of the intermediate amplitudes are 
correlated, e.g., u> and v — uj in Fig. [JJd), there is only one independent probe field 
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frequency uj to be integrated over at fixed detected frequency v. Denoting the signal 
frequency v in diagrams (b2) and (c2) of Fig. finally we obtain the following result 
for the inelastic interference spectrum: 



where equations (jlOap . (|10b[) . and (|10c|) correspond to diagrams (b2), (c2), and (d) in 
Fig.[7J 



4. Conclusion 

We presented a diagrammatic derivation of the double scattering background and 
interference spectral distributions of coherent backscattering from two two-level atoms. 
Although these distributions were known before jTTJ [T21 E5] , here we re-derived them 
in an intuitive way which, we believe, makes the diagrammatic approach to CBS of 
intense laser light from cold atoms more accessible and attractive to the reader. 

It is not only the diagrams that render our present approach intuitive, but also a 
physical interpretation of the building blocks which is given in this work. In particular, 
we established a connection of the 'elastic' building blocks with the response functions 
used to calculate the weak probe absorption spectra by a laser-driven atom [15] . 
Furthermore, the 'inelastic' building blocks include the Mollow triplet, a modification 
thereof proportional to the intensity of the weak probe field, and inelastic response 
functions which are linear in the probe fields, and which emerge due to mixing of the 
laser and probe waves. 

Among the single-atom spectral response functions relevant to CBS, only the 
resonance fluorescence and probe absorption spectra were observed experimentally 
[22l [23] . It remains to be seen whether the other functions can also be measured 
directly in experiments with single quantum scatterers interacting with one strong and 
one weak laser field. Perhaps, strongly driven quantum dots [16] are most suitable 
candidates for such observations. 

The possibility of obtaining the double scattering signal diagrammatically from 
single-atom building blocks suggests two directions of future research, in the field 
of CBS of intense laser light from cold atoms. First, the diagrammatic techniques 
presented in this work will be further developed to assess higher scattering orders of 
CBS. Second, a generalization of this method to realistic atoms having degenerate 
dipole transitions will be another useful application. 
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Appendix A. Expressions for building blocks 

Here we present expressions for the elementary elastic and inelastic building blocks 
(Figs. [4] and [5j respectively). A detailed derivation thereof can be found, for example, 



+ (a+) (0) ( ( T-(^) ( - ) P (+) (0,^ 



(10a) 
(10b) 




(10c) 
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Appendix A.l. Elastic blocks 

The expressions associated with the building blocks of Fig. 0] follow from the 
perturbative (in the probe field amplitude) solutions of the optical Bloch equation 
©. They are given by the first or second entries of the following vectors: 



(<7)(°) = GL, (A.la) 

(a(u)) w = GMAM(<7)(°\ (A.lb) 

= G(-iij)A^(a}^, (A.lc) 

(a(u>))W = GAW(?H>(-) + GA(-'(?(w))W, (A.ld) 
where G(x) is the free propagator given by 

GO) = (x~M)-\ (A.2) 



where M is the Bloch matrix (see equation ([3J)), G = G(0), and the matrices A^ 
describe the coupling of the atomic dipole to the weak probe fields (see equation (|4|). 
For the Bloch vector (a) = ((ct~), (a + ), (a z )), the zeroth-order elastic scattering 
amplitude (er - )(°) is given by (GL)i. Likewise, using (|A. 1|) . one finds all of the 
amplitudes of Fig. SJ 

Appendix A.2. Inelastic blocks 

The expressions for the four elementary inelastic building blocks of Fig. [S] can be found 
by solving the equations of motions for the two-time correlation vectors: 

9i(h,t 2 ) = (a(ti)a-(t 2 )) - (a(ti))(a-(t 2 )), h > t 2 , (A.3a) 
Mtuh) = (a+it^aih)) - (a + (h))(a(t 2 )), h < t 2 . (A.3b) 

By virtue of the quantum regression theorem the vectors gi(ti,t 2 ) of equation (|A.3| 
satisfy equations of motion similar to the OBE @, which can be solved by Laplace 
transformation. The resulting expressions for the inelastic building blocks are given 
as a sum of two terms corresponding to equation (|A.3r ) and (|A.3b ). respectively: 



PW( W ) = [gri-iuj)]! + [g^{iuj)] 2l (A.4a) 

P(-J(w,i/) - ^(w.iw-ii/Jli + ^^w.ti/Jla, (A.4b) 

pM(u,u) = + ^ +) (u,w-iu)h, (A.4c) 

P^{u,v) = $P(u,-w)]x + [f?\w,iv)h, (A.4d) 
where the Laplace transform solutions read: 

l, 2 (z)^=G(z)^(0), (A.5a) 

Hij(u,z)M =G(z + tuj)[AM^l(z)+g[+\u,0)}, (A.5b) 

| 1)a (u,«)(-) =G(z-iu)[A(-^l(z) + g^(u,0)}, (A.5c) 
fi >2 (a;,*)< 2 > = G( 2 )[A(+>g( W)Z ) + A<->|g(w,*) 

+ ^3(w,0)]. (A.5d) 
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The initial conditions in the right hand side of equation (|A.5[) result from setting 
h = t 2 in CO}: 

S< 0) (0) = -i A<+> <<?> W + L x - (a) C°> (a~ > (0) , ( A.6a) 

$\0)=iA^(*)W+L 2 - (a)^(a+)^, (A.6b) 

5^(^,0) = -iA< + ><*M> (±) - (<?K» (± V> (0) 

-<*> ( °V~H> (±) > (A.6c) 
#>(«,<)) =iA(")(a(u;))( ± ) - (a(u,)> (± V> (0) 

- <^> (0 V + H) (±) , (A.6d) 
sf^O) = ^A(+)<aH> (2) - <aH> (2) (O (0) 

- <^> (0) <^M) (2) - (a(w)> (+ V~(«)) (_) 
-<<?M) ( -VM> (+ \ (A.6e) 

-(a)(°)(a+M)( 2 )-(a(a;))(-)(a+M)W 
-(?H>W(^H}'-' (A.6f) 

with Li = (0,1/2,0) T , and L 2 = (1/2,0,0) T . 

Appendix B. Relation of the building blocks to the results of Mollow 

Apart from notation, the elastic building blocks are equivalent to the results obtained 
by Mollow, in connection with the study of stimulated emission and absorption of 
a weak probe field by a coherently pumped two-level atom [15] • Since Mollow's 
derivation and definitions differ from ours, we here establish the relation between 
his and our results. Upon the replacements 

k — > 27, k' 7, Alu A, 

Av -> lu, 2|A£ | 0, A'£„ -)■ u, 

where the left and right hand sides of each substitution correspond to the results of 
|15j and of our present work, respectively, we obtain the following identities: 



a 


= <O (0) , 


(a)* 




(B.l) 


5a+ 


= v(o--(u))(-\ 






(B.2) 


5a,- 


EE«(a+H)(-), 


(8a-)* 


= «>"(«)>« 


(B.3) 


Scto 


= Av 2 (a-{u)) {2 \ 


(Sao)* 


= 4( W *)V + M) (2) - 


(B.4) 



We note that the right hand side of each identity is equal, up to a pref actor, to 
the elementary building block from Fig. |4] For completeness, we also provide the 
expressions on which the elementary building blocks depend implicitly: 



n-m= (a z ) {0 \ rj = 2v(a z (uj)) ( -\ (B.5) 

Sn = 2\v\ 2 (a z (cu)) (2 \ (B.6) 
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where, as in equations (IB.ip ~ dB.4p . the left and right hand sides of each identity refer 
to our present results and those of [TS], respectively. 

As shown in [TS], the rate W' of absorption/stimulated emission of the probe field 
is given by 

W = -i\'£'*8a+ + i\' £' Q (8a + )* . (B.7) 

Hence, the identities (|B.2j) imply our interpretation of the building blocks in Fig.|4](c,d) 
as describing a stimulated emission of the probe field in presence of the pump field 
(whether the field is actually emitted or absorbed is defined by the sign of W' [15 ). 
Analogously, the rate <5W of absorption/stimulated emission of the pump field reads 

8W = -iX*££5a + iX£ {5a )* . (B.8) 

Therefore, due to the identities (|B.4j) . we interpret the amplitudes (a^ (oj))^ (see 
Fig- El(g,h)) as describing stimulated emission of the pump field in presence of the 
probe field. 

The physical interpretation of the function 5a- (((5a_)*) was not considered by 
Mollow. But, by analogy with the above results, and using the fact it represents a 
correction to the negative- (positive-)frequency scattering amplitude for a pumped 
atom probed by a positive- (negative-) frequency field, we deduce that (5a_ ((fo_)*) 
and, consequently, (er+ ({a~ see Fig. IUe,f), describes a phase- 

conjugation due to non-linear mixing of the pump and probe fields. 
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